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Abstract 

In order for surface scattering models to be accurate they must necessarily satisfy energy conserva¬ 
tion and reciprocity principles. Roughness scattering models based on Kirchoff’s approximation or 
perturbation theory do not satisfy these criteria in all frequency ranges. Here we present a surface 
scattering model based on analysis of scattering from a layer of particles on top of a substrate 
in the dipole approximation which satisfies both energy conservation and reciprocity and is thus 
accurate in all frequency ranges. The model takes into account the absorption in the substrate 
induced by the particles but does not take into account the near-held interactions between the 
particles. 
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1. INTRODUCTION 

In this paper, we report a model for scattering at an interface separating two homogeneous media. 
The aim is to derive a BSDF (Bidirectional Scattering Distribution Function) that accounts for 
both specular and diffuse components of light scattered both in transmission and in rehection, and 
also absorption. The initial motivation of this work is to develop a model that can be used to 
analyze multiple scattering in a photovoltaic cell. Such a cell is a multi-layer system comprising 
an absorbing active layer (e.g. silicon or GiGS) deposited on different materials that can be used 
as electrical contacts, back-rehectors, antirehection coatings, etc. Most active materials used for 
photovoltaic cells are semiconductors whose refractive index takes large values so that accounting 
for total internal rehection is very important. A basic mechanism that is often used to increase 
the absorption is to scatter light in order to trap the light in the absorbing medium. Most designs 
are performed either using a periodic model for the roughness or by trial and error using randomly 
rough surfaces. 

In this paper, we explore the issues raised by the modeling of scattering and absorption in such 
a structure using a radiative transfer equation approach. For a multi-layer geometry, an adding 
doubling approach to the solution of the radiative transfer equation is appropriate. Hence, the 
basic tool needed, is a matrix accounting for the BSDF of the interface. The key issue as far as 
photovoltaic applications are concerned is that the BSDF needs to fulhl energy conservation with 
an accuracy better than 1%. Indeed, the whole point of the modeling of a photovoltaic cell is to 
gain a few percent in the absorption. This sets the standards required for the accuracy of the 
model. Not less important is the requirement of reciprocity. Indeed, when a cell is designed in 
order to scatter light to couple incident propagating light to trapped light in the cell, the same 
roughness can couple back trapped light to propagating modes with the same efhciency owing to 
reciprocity. Hence, if light can be coupled to guided modes, then guided modes can be coupled 
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back to propagating modes and escape as shown in Fig. Thus, optimizing the right coupling 
between guided modes and propagating modes is far from trivial and requires an accurate model 
that accounts properly for reciprocity. 

Currently available models [DEli for scattering by rough surfaces are focussed on modeling 
properly the angular scattering pattern and describing properly the diffuse and collimated reflection 
factors. Satisfying reciprocity and energy conservation is not the major issue for most of them. As 
a matter of fact, models based on Kirchhoff approximation, phase perturbation, etc do not satisfy 
energy conservation and reciprocity with good accuracy over all frequency ranges. For energy 
applications, what is needed is not just an exact angular description of the scattering but also 
a correct modeling of the balance between collimated and diffuse scattering as well as a correct 
balance between scattering and absorption. When dealing with light propagating in scattering 
media such as gases, particles or paper for example, the isotropic approximation for the phase 
function is very often used @]. Here, we establish a similar approximation while accounting for 
the presence of interfaces. Inasmuch as the system is in a multiple scattering regime, the final 
result does not depend much on the details of the phase function. Instead, accounting properly 
for absorption and scattering and satisfying energy conservation is critical. Hence, approximations 
such as Milne-Eddington, diffusion approximation, [SIS] etc are extremely useful. This discussion 
sets the landscape: we seek a model that allows introducing scattering and absorption with enough 
parameters to control the balance between collimated and diffuse, and between scattering and 
absorption. As we are looking for a multiple scattering regime, the accuracy in modeling the exact 
angular behavior is not very important. By contrast, it is of critical importance that the model 
satisfies energy conservation and reciprocity. 

In order to establish such a model, we start from the scattering of a single dipolar scatterer located 
at a certain distance from an interface. For such a scatterer characterized by its polarizability 
a such that its dipole moment is proportional to a and the incident field on the scatterer, we 
can derive explicitly the scattered field in an electrodynamics framework. This solution satisfies 
energy conservation and reciprocity. It does also account for the coupling of light into guided 
modes in the medium with high refractive index. Finally, it also accounts for absorption induced 
in lossy substrates by the near-held produced by the scatterer. This elementary result, which was 
extensively discussed in Ref. [7] will serve as a building block for our model. In this paper, we take 
a further step in order to develop a model for the scattering. We need to establish a model for the 
coherent modihcation of the rehection and transmission factor for the collimated beams as well 
as account for diffuse scattering and absorption. The available model parameters are as follows: 
the real and imaginary part of the polarizability, the number of scatterers per unit surface, and 
the distance between the scatterer and the interface. In the low density regime, our model will be 
expected to ht with an exact numerical solution of the problem of scattering by N scatterers. For 
larger density, near-held and correlation effects are expected to become important. Our goal in 
this paper is not to capture these effects accurately. Instead, we only look for a model that allows 
introducing a degree of scattering while preserving energy conservation and reciprocity. 
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Figure 1: Light trapped in a layer can undergo both coherent reflection and transmission as well as multiple scattering 


The paper is arranged as follows: In Sec. a brief description of scattering from a particle in a 
homogeneous medium is given which introduces the terminology used in this paper. This is then 
extended to the case of scattering from a layer of particles in a homogeneous medium in Sec. 
where we use the effective index method based on Maxwell Garnet theory to analyze the coherent 
scattering from the particles. We then use this theory in Sec. to obtain the BSDF of a surface 
with a layer of particles on top of it serving as source for the scattering function and compare it 
in Sec. [^with exact numerical simulation of scattering from a layer of particles on top of a finite 


surface. Finally, in [Appendix A| we show the microscopic derivation of coherent scattering from 
a layer of particles in a homogeneous medium where we take into account the multiple scattering 
between the particles using the mean-field theory but ignore near-held interactions and correlations 
(and hence recurrent scattering) between them, thereby establishing the equivalence between the 
effective-index model and coherent scattering in the mean-held approximation. 


2. Scattering from a single particle in a homogeneous medium 
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Figure 2: Configuration for finding the scattered field at any point (r, z) from a particle with center located at 
(0, 20) in a homogeneous medium of dielectric function ei. 


In the conhguration shown in Fig. consider a scatterer located at (0, zq) in a homogeneous 
medium of dielectric permittivity Si. For simplicity of analytical description we assume that the 
surrounding medium is non-absorptive so that ei is real. The incident held at position (r, z) 
is chosen to be a planar wave with amplitude unity at the position of the particle, and with a 
wavelength much greater than the radius of the particle a. The incident held can be written as: 

E,(r, z) = e, _ ( 1 ) 

where is a real vector denoting the component of the incident wave vector parallel to the xy 
plane; k^i is related to kjj^^ as: 

kj = ei(a;/c)2 = kj^ + (2) 

where uj is the frequency of incident radiation and c is the velocity of light, is the polarization 
of the incident wave: 

[s = (Ox + (—l)y-1-Oz) s/TE— polarization 

Bi = < II (3) 

I^Pi- = (fczix -I- Oy -I- kjJ^j, z)/fci p/TM — polarization 

For a TM polarized wave travelling in the positive z direction we would have = pi+ where 
Pi+ = (-fczix -k Oy + z)/ki. 

The incident held polarizes the particle such that its dipole moment is given by [8]: 

Po = eoaoEexc(0, zo) . (4) 


where, Eq is the electric permittivity of free space, is the polarizability of the spherical particle 
given by: ao = 47ra^(ep — ei)/(ep -I- 2 ei); Sp is the dielectric permittivity of the particle, and Eexc 
is the exciting held external to the particle and comprises of not just the incident held E^ but 
also the radiative reaction held (scattered held) [8] from the particle. If we use the (dyadic) Green 
function which relates an electric-dipole source p(r') at a position r' to the electric held E(r) at a 

position r through the relation E(r) = G(r,r')p(r') we can write the exciting held Eexc as: 

Eexc(r) = Ei(r) -f Go(r, r')po(r') 


o 


( 5 ) 







where in the limit of a{uj/c) —)■ 0 the Green’s function Go(r,r') is given by [5]: 


Go(r,r') « i 
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( 6 ) 


with I being the unit dyad. If we express the dipole moment po in terms of the incident field Ei 
as: 

Po = £oQ!effEj(0, Zo) • (7) 


from Eq. mu and 1^ we get: 
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1 — /(67r)ao 


( 8 ) 


It must be noted that the expression for aeff in Eq. [^is valid in the limit aw/c —)■ 0. For larger 
particles an expression for apff can be derived from Mie theory of scattering from a sphere and is 

given by mn: 

Gtt 


^eff — 


kl^{CE-i) 


(9) 


where the coefficient Cb is: 
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( 10 ) 

with gj, and gm being non-dimensional factors given by y/Ii{ujlc)a and yC^{ui/c)a respectively. 
The scattered field from the polarized particle with dipole moment po given from Eq. can be 
written in Sipe’s formalism [12] as: 


E(r,z)=/ ^F(k||,z-zo)e*'^" 

Jkii 


( 11 ) 


where, 
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(13) 


(14) 


The key point here is that this form of the scattered field satisfies energy conservation and reci¬ 
procity. We now proceed to derive a model for interface scattering using this result as a building 
block. 


3. Effective index model 


For a layer of particles located in a homogeneous medium it is possible to derive the specular 
reflection and transmission coefficients of an incident planar wave by finding the mean scattered 


field from the particles (as shown in Appendix A). However, while this derivation throws light 
on the coherent effects in scattering by the particles, it is not convenient to adopt this approach 
for finding the coefficients when interfaces are present in the surrounding medium. To account 
for multiple reflections between the interfaces and the layer of particles, we develop an alternative 
approach based on the effective index theory to derive the reflection and transmission coefficients 
for the case of a layer of particles in a homogeneous medium, show that the expressions for the 
coefficients are equivalent to that derived from the mean-field theory, and look for ways to extend 
it to the case where interfaces are present. 


















In this method we use the Maxwell Garnet theory to replace the layer of particles by a thin film 
of arbitrarily small thickness and with a refractive index such that the total polarization in the 
film is the same as that of the layer of particles. A plane wave would thus interact with such 
a film in the same way as it would with the layer of particles. If eeS,x, £eS,y and eeB,z are the 
dielectric properties of the film along the x, y and z directions respectively, E^, Ey and Ez are 
the components of the electric field above the effective index layer in the incident plane (here the 
boundary conditions necessitates that Ex and Ey are the same in the film and outside), p is the 
density of particles (units of m“^) in the layer and deff is the thickness of the film, equating the 
total polarization per unit area of the film to the total polarization in the layer of the particles we 
get: 
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This gives us the dielectric properties of the film as: 


^e.B^x — S'eff^y — T P^&bI deff) 

^eB,z — 1/ (f P^eB/{^ideB^^ • 


(15) 


The polarizability of the particles, Oeff, is isotropic for the case of particles located in a homogeneous 
medium (this condition, however, gets relaxed when interfaces are present as discussed in Sec. |^. 
The reflection and transmission coefficients of a plane wave incident on the thin film is given by: 


and 


reS = 


^eB 


j' _ rp 

(16) 
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where keB 
given by: 


y^£eff,s/(w/c)2 


r is the Fresnel reflection coefficient which for a TE mode is 
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and and are the Fresnel transmission coefficients given by (for TE mode): 
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we can show that: 
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and 
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( 22 ) 


The corresponding coefficients for the TM mode are slightly more complicated on account of 
the different dielectric properties along the x and z directions, and hence the coefficients do not 
reduce to simple expressions like in Eq. and The form of the Eresnel coefficients to be 
used in Eq. and for the TM modes are given in [Appendix B| The expressions for the 
reflection and transmission coefficients shown in Eq. are equivalent to the expressions for 


















the coefRcients that are derived by considering the mean scattered field of individual particles in 
the layer (see [Appendix A I thereby implying that the coherent scattering by a layer of particles 
can be equivalently modeled by a thin film of a fictitious material whose dielectric property is 
given by Eq. |15| Thus from the energy conservation statement across a thin dielectric film it is 
possible to arrive at an expression for energy conservation across the layer of particles, as follows. 
The energy conservation statement across a thin dielectric film reads: 


Re(Si • z) = Re(S 2 • z) + / Re 


IUJ£q , <->■ 


(Eeff- l)Ed(r,2;).E2(r,z) 


dz 


(23) 


where Si and S 2 are the Poynting vectors of the coherent field at interfaces above and below 
the layer of particles in the homogeneous medium respectively and the integral term in Eq. |23| 
represents the extinction in the coherent field across the layer of the particles, with Ed(r, z) being 
the field inside the thin-film, is a diagonal tensor with elements (£eg,x, £eS,y, £eS,z)\ indicates 
the complex conjugate and dz being a unit element of the effective index layer along the z axis. 
The integral is over the thickness of the thin dielectric film. Eq. can be shown to reduce to the 
form: 
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(keffP + IteffP - 1) = / Re 


IUJSq 


(eeff-l)Ed(r,z).E*(r,z) 


dz 


where, reS and teS are the reflection and transmission coefficients from Eq. |21| and [^ 
the values from Eq. 15 and taking the limit dcff/ (poio) —>■ 0 we can write: 


(24) 

Substituting 
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(25) 


where Eiii(r, 2 : 0 ) is the field outside the film thus comprising of both the incident and reflected 
fields. 

The right hand side of Eq. [^is the extinction in the coherent component of the incident radiation, 
which comprises of both the diffusely scattered radiation as well as the flux absorbed by the particles 
in the layer. This can be seen by comparing the expressions for the power scattered and absorbed 
by a single particle. The power absorbed by a particle in a homogeneous medium is given by [5]: 


fEabs = 


weo 
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The field scattered by a dipole in a homogeneous medium is given by Eq. E] from which we can 
show the power scattered to be m- 
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(27) 


From Eq. [M 
as: 


25l and we get the energy conservation statement across the layer of particles 
kiz 


2u;/io 


(|reff| + |/eff| - 1) = /0R4ca + pWabs 


(28) 


with pWsca, accounting for the total diffuse scattered power from the particles. It is worth pointing 
out explicitly that the illuminating field Ein includes not just the incident field Ei from Eq. 
but also the coherently scattered field from the particles. This can also be seen in the microscopic 
description of scattering from the layer of particles shown in |Appendix A| where the mean scattered 
field from the particles has been shown to account for the coherently reflected field. The dipole 
moment of the particles in the layer, po, can thus be written as: 


Po = eoQ!o(E exc + Erefl) 


(29) 


where Ej-efl is the coherently reflected field from the layer of particles and Eexc is given from Eq. 

As such, while the diffuse scattered power pWsca. in Eq. seems to resemble the independent¬ 
scattering approximation El: we point out that the model does take into account the interaction 
between the particles in the mean-field approximation. However the model does not take into 
account the near-held interactions and correlations in the position of the particles. The model 
inherently satishes reciprocity principle since it makes use of an analytic solution of the electro¬ 
magnetic scattered Held from a particle which is known to satisfy reciprocity [5]. This will be 
shown explicitly in Sec. |4.2| 






















4. Analysis of coherent scattering from a mono-layer of particles above a substrate 


As explained in Sec. [l]the intention of this work is to develop a BSDF of a surface which satisfies 
energy conservation and reciprocity. With this aim in mind and with the results from the previous 
sections we now attempt to find the BSDF of a surface with a layer of particles on top of it acting 
as the source of scattering as shown in Fig. ia). Scattering from a single particle on top of a 
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Figure 3: (a) Configuration for finding the scattered field from a layer of particles located at a distance zq in a 
homogeneous medium of dielectric function ei on top of a substrate with dielectric function £ 2 . Surfaces 1 and 2 
are located just above the position of the particles at 2 = Zq and just below the surface of the substrate at 2 = 0“ 
respectively, (b) In order to take into account multiple reflections between the layer of particles and the substrate 
the layer of particles is modeled as a thin film of arbitrarily small thickness with effective dielectric function Egff 
and located at a distance 20 from the substrate 


substrate has been discussed in detail in Ref. [7] . The main results from this analysis is presented 
in Sec. 4.1 and the results are utilized to extend the analysis to the case of a layer of particles on 
top of a substrate in Sec. |4.2| 


4-.1. A single particle above a substrate 


Consider the case of a single particle located at (0, zq) above a substrate comprising the half-space 
z < 0. In the vicinity of the substrate the radiation reaction field from the particle discussed 
in Sec. [2] will be altered due to reflection of the scattered field from the substrate below. The 
effective polarizability from Eq. [^will thus have to be modified to take into account this reflected 
component. The procedure to arrive at this modified form of Oeff is similar to the one detailed in 
Sec. with one difference - the effective polarizability will now be a tensor. The polarizability 
tensor Oeff is diagonal with elements (a^, a||) and with ^ a^g as the presence of 

the surface breaks the rotational symmetry. The exciting field Egxc in Eq. [^is now given by: 

Eexc(r) = Ei(r) -k [Go(r, Tq) -k Gr(r, ro)]po(ro) (30) 

where G^ is the (dyadic) Green’s function which accounts for the reflection component of the 
scattered field from the interface and can be calculated by the Fourier transform [7]: 

Gr(r,ro)=/ -^Gr(kll;z-zo)e*’" ■ (31) 

7_oo (27ry 
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(32) 


Here, r™ and r™ are the Fresnel reflection coefficients for the TE and TM modes respectively at 
the interface between the homogeneous medium of dielectric function £1 and the substrate. The 
dipole moment of the particle is given by: 


and from Eqs. [6 
form: 
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and 


33 


Po = £oQ;effEi(0,zo) . (33) 

it is possible to show that the effective polarizability Seff takes the 
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where I is the unit dyad. 

The scattered power Wgca from a particle with dipole moment po given in Eq. [33| is no longer 

(R) ' -' 

independent of the orientation of the dipole. If Wsca denotes the scattered power from the particle 
at a surface just above the position of the particle at at z = 0+, and if wira denotes the scattered 
power from the particle at a surface just below the surface of the substrate at z = 0“, expressions 
to calculate Wsca"^ and 114^'^ as an integral over k|| have been derived in Ref. [7]. These are given 
by: 
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(38) 


Here and t™ are the Fresnel transmission coefficients for the TE and TM modes at the 
interface between the homogeneous medium of dielectric function Si and the substrate respectively. 
The total power scattered from the particle lEsca can be computed Wsca = W'sca'^ + Hsca ■ Note 


that Eq. 37 includes interference between the light scattered upward from the particle and the 
light scattered downward from the particle and subsequently reflected from the substrate below. 
Another important point to be noted here is that the proximity of the particle to the substrate 
induces contributions from large wave-vectors (fe|| > w/c) to the scattering power M4 
will be seen as absorptive losses in the substrate [7]. 


which 


4--2. Layer of particles above a substrate 

With the results from Sec. [^and Sec. |4.1| we are now in a position to describe scattering from a 
layer of particles above a substrate. To obtain the coherent reflection and transmission coefficients 
for the configuration shown in Fig. [^a) of a plane wave incident on a substrate with particles on 
top of it, we adopt the procedure detailed in Sec. and replace the layer of particles with a thin 
film of thickness deff as shown in Fig. [^b). To make sure that the total polarization is the same in 
both the cases and noting the anisotropic nature of the polarizability of the particles, the dielectric 
properties of the thin-film are now given by: 


^efF,x — ^ef^,y — (f “b P^efil ^eff) 

£es,z = 1/ (1 - pale/i^ides)) ■ (39) 

With these changes and employing standard recursion relations m or transfer-matrix method 
|16j we can compute the coherent reflection and transmission coefficients for the multi-layered 
configuration shown in Fig. ib). The energy conservation statement from Eq. is still valid 
for this configuration, with the dipole moment of the particles po to compute Wgca being given by 
Po = where the illuminating field Ein is taken to be the field above the thin-film in Fig. 

Sb). This field accounts for the multiple reflections between the particles in the layer as well as 
with the substrate below. 


In radiometry the BSDF is frequently used to characterize the scatter of optical radiation from a 
surface as a function of the angle of the incident and the scattered beam. It is given by the ratio 
















of the scattered radiance to the incident irradiance and has units of sr“^. The terms bi-directional 
reflectance distribution function (BRDF) and bi-directional transmittance distribution function 
(BTDF) are used when referred specifically to the reflected and transmitted scatter respectively. 
If (j)) and (/>) are the specific intensities of the scattered beam (units of Wm“^sr“^) at 
interfaces 1 and 2 respectively shown in Fig. |^a) where (j) is the azimuthal angle and fi = cos 9 
with 9 representing the polar angle of the scattered beam, we can express the BRDF and BTDF 
for the surface shown in Fig. j^a) on which a coherent beam given in Eq. [^is incident upon it as: 


BRDF(/x, (j)) 
BTDF(/x, (/>) 


{2u;fj.o)I+{p,<j3) 
Re[(E,x (k* xE*)).z] 

Re[(E,x (k* xE*)).z] 


(40) 


To arrive at an expression for and /“(/r, </)) in Eq. 40 in terms of the electromagnetic 

description given in the previous sections we compare the expressions for the reflected flux above 
the surface. If Flux*-^^ denotes the radiative flux at surface I in Fig. [^a) it is given in terms of 
as: 


Flux(^) = 




(41) 


^=0 J ^—0 


But from Eq. [28] and the description given in Sec. 4.1 this must be equivalent to the diffuse 
scattering from the layer of particles at surface 1 which is given by: 


Flux(^) = 


(42) 


with the integral over fc|| in Eq. 35 extending from 0 to Y^(a;/c). Utilizing the cylindrical 


transformation (i^k|| = dk\\ dcj) in Eq. 35 and noting that fc|| = ki sin0 we get d^k|| = 


Thus Eq. 35 can be written as: 

P‘27t p 1 

^sS{z = z^)= / M+{fi,(l))kfndfj.d(j) 

Jo Jo 


(43) 


with 




(2,)> I 


F,,(k||;zo) X (k* xFF(k||;zo)) 


.z 


From Eqs. \4^ and [43l we get: 




Substituting this in Eq. 40 we get for the BRDF: 
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Similarly, the expression for BTDF reads: 
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The reciprocity relation for this model can be confirmed by comparing the differential scattering 
cross section in the scattering direction {9sc,4>sc) for an incident beam at angle (0^, (j)i) in the config¬ 
uration shown Fig. with that obtained by interchanging the incident and scattering directions. 
That is, for reciprocity to hold true we need m- 
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where represents the differential scattering cross section. In terms of the BRDF, Eq. 
ail 

equivalent to: 


BRDF i (/Xgc, 4>sc) Msc — BRDF sc ■, f^i 


49 


is 


(50) 


where BRDFi(^sc 5 ^^'sc) [BRDFsc(/ii, represents the BRDF from Eq. 46 for light incident in 
the direction (0^, (j)i) [(0scj </^’sc)]i Mac = cos^gc, and /i^ = cos0i. A similar expression can also 
be written for the BTDF. We confirm the reciprocity relation for our model by taking two pairs 
of arbitrary angles for the incident and scattering directions as shown in Table and finding the 

+ ^(6'sc,()';6'i,<)>) fora 


error calculated as 


fixed azimuthal angle ^ = 0. The values shown are for the average of TE and TM polarizations. 
The material parameters in the configuration are arbitrarily chosen with the polarizability of the 
particles ao = 6.06 x 10® nm®, dielectric function of the substrate £2 = 3.91 + 1.2f, incident beam 
of wavelength A = 300 nm and filling fraction of the particles / = 0.15. The non-dimensional 
quantity / is related to the density of the particles p as / = pira^. The reciprocity relation is 
observed to be satisfied up to floating point precision. 


9i, 9sc 

error 

40°, 74° 


58.3°, 85.1° 

T 

0 

1—1 


Table 1: Error in satisfying the reciprocity relation in our model for two pairs of arbitrarily chosen angles 


Noting that the integral over the azimuthal angle (f) in Eq. can be analytically obtained, the 
scattered flux in Eq. |^can be alternatively expressed as: 

eW^) = f (51) 

J fl 

we have ~ fo^ (j))kld(j). We can use this to define the 


where, from Eq. 


41 


and 


BRDE alternatively as: 
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BRDF(p) = 


(2wpo) in, (t))kld(t) 


(52) 


Re[(E, X (k* X E*)).z] 

We now have all the tools required to describe the utility of the scattering model developed so far. 



Figure 4: BRDF from Eq. |52| for the configuration shown in Fig. [^with the following parameters: £2 = 3.91 + 1.2i; 
A = 300 nm; / = 0.15; ao = 6.06 X 10® nm® and (a) zq = 85 nm, (b) zq = 100 nm, and (c) zq = 135 nm. The 
parameters for Fig. (d) and (e) are the same as in (b) and (c) respectively except for ao which is increased to 
9.87 X 10® nm® 


As a demonstration, consider the particular case of scattering from the surface shown in Eig. j^a) 
with the material properties taken to be the same as that chosen for showing reciprocity relation 
in Tableand assuming normal incidence. The BRDF from Eq. 52 is plotted in Eig. |^a)-(c) as a 
function of the scattering angle 9 for varying position of the particles above the substrate, Zg- What 
we observe is that by increasing zg scattering lobes with large angle scattering can be obtained. 
Likewise, keeping zg constant if we increase the polarizability ag of the particles (or equivalently 






























the filling fraction /) we can increase the scattering cross section of the particles. This can be 
seen by comparing Fig. |^b) and |^d) (or similarly Fig. Qc) and|^e)) where zq has been kept 
constant but ao has been increased. We observe that the magnitude of scattering has increased 
but the shape of the scattering profile has been retained. Thus by varying Zq and ao (or /) we can 
change both the scattering angle as well as magnitude of scattering and thus obtain a wide range of 
scattering profiles. This can be used, for example, to analyze the effect of scattering in increasing 
the absorption of light in the active layer of a multilayered solar cell. By tuning the parameters 
Zq and ao a wide range of scattering profiles can be explored so as to optimize absorption in the 
active layer. 


5. Comparison with exact electromagnetic simulations 


The scattering model developed in this work can be compared with exact electromagnetic simula¬ 
tion of scattering by a random layer of spherical particles located on top of a finite substrate and 
averaged over many realizations. Details on the procedure for the exact electromagnetic simula¬ 
tions can be found in Ref. [10]. It must be noted that the simulation will include the near-field 
interactions between the dipoles and hence is expected to diverge from the scattering model devel¬ 
oped in this work for large filling fractions of the particles. In particular, it has been shown in an 
earlier work [18| that near-field interactions and correlations in positions of the particles become 
important when the filling fraction of the particles is greater than approximately 5 %. 

Here we show a comparison between (i) specular reflection as predicted by the effective index model 
in Eq. and (ii) diffuse scattered power from Eq. [^(both quantities normalized with the incident 
power) with the corresponding quantities observed from exact electromagnetic simulations. This 
has been shown in Table where we consider a particular case of scattering by a layer of particles 
lying on top of the substrate (zq = a) with arbitrarily chosen parameters: ao = 2.88 x 10® nm®, 
S 2 = 3.91 -I- 1.2i, A = 300 nm, and for different filling fractions 5% and 15%. As expected the 
results deviate significantly for the higher filling fraction. The deviation for higher / should not be 
viewed as a setback for our model since it was never our intention to model the multi-scattering 
problem accurately but rather we wished to formulate a scattering model which gives importance 
to satisfying energy conservation and reciprocity principles. 


a) 

filling fraction = 5% 

o 

O 

II 

O 

o 

CO 

II 

6»i = 60° 

ao = 2.88 X 10® nm® 

Exact 

Approx 

Error(%) 

Exact 

Approx 

Error (%) 

Exact 

Approx 

Error(%) 


SR without particles 

0.1199 

0.1199 

- 

0.1564 

0.1564 

- 

0.3350 

0.3350 

- 

TE 

SR with particles 

0.0290 

0.0313 

7.9 

0.0420 

0.0480 

14.3 

0.1600 

0.1805 

12.5 


diffuse scattered power 

0.1160 

0.1261 

8.6 

0.1240 

0.1349 

8.8 

0.1310 

0.1388 

7.3 


SR without particles 

0.1199 

0.1199 

- 

0.0872 

0.0872 

- 

0.0062 

0.0062 

- 

TM 

SR with particles 

0.0290 

0.0313 

7.9 

0.0280 

0.0317 

13.2 

0.0120 

0.0115 

4.2 


diffuse scattered power 

0.1160 

0.1258 

8.4 

0.1000 

0.1112 

11.2 

0.0890 

0.1049 

17.9 


b) 

filling fraction = 15% 

o 

O 

II 

Oi = 30° 

6»i = 60° 

ao = 2.88 X 10® nm® 

Exact 

Approx 

Error(%) 

Exact 

Approx 

Error(%) 

Exact 

Approx 

Error(%) 


SR without particles 

0.1199 

0.1199 

- 

0.1564 

0.1564 

- 

0.3350 

0.3350 

- 

TE 

SR with particles 

0.0160 

0.0094 

41.3 

0.0130 

0.0226 

73.8 

0.0800 

0.1533 

91.6 


diffuse scattered power 

0.1750 

0.2476 

41.7 

0.1890 

0.2577 

36.4 

0.2020 

0.2503 

23.9 


SR without particles 

0.1199 

0.1199 

- 

0.0872 

0.0872 

- 

0.0062 

0.0062 

- 

TM 

SR with particles 

0.0148 

0.0094 

36.5 

0.0025 

0.0052 

108 

0.0200 

0.0249 

24.5 


diffuse scattered power 

0.1800 

0.2477 

37.6 

0.1790 

0.2705 

50.2 

0.1900 

0.3182 

67.4 


Table 2: Normalized specularly reflected (SR) power from the effective index theory, and the diffuse scattered power 
from Eq. |51| (marked ‘Approx’) are compared with values from exact electromagnetic simulations (marked ‘Exact’) 
for the configuration shown in Fig. I^a) for different incident angles (Oi = 0°, 30° and 60°) and for two different 
filling fractions (a) 5 % and (b) 15 %. The error between the simulated and approximate values are expressed in 
terms of percentage of the exact values. Specularly reflected power from the substrate in the absence of the layer of 
particles is also presented for reference. 






























6. Conclusion 


With the motive of exploring different scattering profiles to be used in optimizing absorption 
in the active layer of a multilayered solar cell, we have developed a scattering model based on 
scattering from a layer of particles on top of a substrate which satishes both energy conservation 
and reciprocity and takes into account the particle-induced absorption in the substrate. In course 
of developing such a model we have shown that the effective index model based on Maxwell Garnet 
theory can sufficiently describe the coherent scattering from the layer of particles. While the intent 
of developing this model is not to accurately describe the scattering from a layer of particles, we 
have shown that in the regime where near-held interactions and correlation in the position of the 
particles are negligible the model sufficiently approximates both the coherent scattering as well as 
the diffuse scattering from a surface with a layer of particles on top of it. 
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Appendix 


A. Microscopic analysis of scattering from a mono-layer of particles in a 
homogeneous medium 


Incidentfield 

z - axis 



Figure A.5: Configuration for finding the scattered field at any point (r, z) from a layer of particle with centres 
located at (r„,zo) in a homogeneous medium of dielectric function ei. 


Here we consider a homogeneous medium containing a monolayer of a large number N of identical 


spherical particles with center of the particles at {ri,zo), {r 2 , Zo)...{r]\f, zq) as shown in Fig. A.5 


We assume that the particles are uniformly distributed so that the medium is statistically homoge¬ 
neous. Propagation of an electromagnetic wave in such a medium is characterized by the presence 
of a mean held (E(r, z)) corresponding to a plane wave in the statistically homogeneous medium 
and a huctuating component corresponding to the huctuating dielectric function |14j . Thus when 
an incident held given by Eq. impinges on the layer of particles, the scattered held from the 
particles comprises of both a coherent component as well as a diffuse component (in contrast with 
the case of scattering by a single particle discussed in Sec. [^with only a dihuse component). 

The mean held (E(r, z)) is given by: 


(E(r,z)) = J P(ri,r 2 , ...,T'Ar) E(r, z) (fri tfr 2 ...(frN 


(A.l) 


where, P(ri, r 2 ,..., rjy) is the probability density of hnding particles at positions {ri, zq), (r 2 , Zq) -^o) 

and E(r, z) is the held at position (r, z) given by: 


E(r, z) = Y^ E(”) (r - r„, z - zq) 


(A.2) 


1 n 







where El"^(r — r„, 0 — zq) denotes the scattered field at (r, z) from a particle located at position 
The scattered field from a particle is given by Eq. 12 so that Eq. A.2 can be written as: 

d^kfi 


E(r,z)=^ 


'kf, (27r)2 


F„(k* z-zo)e*'^ii 


f*—"E tj.) 


(A.3) 


The dipole moment of the particle located at (r„, zq), Po"\ can be written as: 

p[,”^ = £oaeffE[j"^ (A.4) 

where e[^^ is the field incident on the particle at position (r„,zo) comprising of not just the 
incident field E^ but also the scattered field at (r„, Zq) from the other particles in the layer. Thus 
the scattered field from Eq. |A.3| reduces to: 

" d^kf, j 


E(r,z)=^ 


L (27r)2 2“‘^®^Hc 




2 gi/c. 


(z-zo) 


kzl 


-el 


k.,.(r-r„ 


SS + P 1 + P 1 + 


•E 


(n) 


ill 


The mean field is thus given by: 

d^kfi 


(E(r,z)) = 




0 


2 gifc^i(z-zo) f 
- ell- 


St 


SS + P1 + P1 + 




(A.5) 


(A.6) 


For brevity we simplify Eq. A.6 for TE polarized waves. Similar arguments can be adopted for 
analyzing TM polarized waves. We first make an assumption that there is no correlation between 
the positions of the particles i.e., we write P(ri,r 2 , ■■■, r^v) = Pi{'i'i)P 2 {f‘ 2 )---PN{f‘N) where Pi{ri) 
is the probability of finding a particle at position {ri,zo). For non-correlation we have Pi(ri) = 1/S, 
(i = 1,2, ...TV) where S is the area in the xy plane where the particles are located. This assumption 
is valid in the dilute regime where the filling fraction of the particles given by / = Ntto^/S « 0 
where N is the total number of particles. For higher filling fractions correlation between the 
positions of the particles have to be taken into account as detailed in Ref. M- We also take 


E|"UiVE,(r,z) 


(A.7) 


where K is an unknown quantity which accounts for the variation of the illuminated field from the 
incident field Ei(r, z) du e to the scattered field. For an incident TE wave, taking the definition of 

and substituting lims^oo ^ = lirng^oo ■^-^d(k||^^- 


A.l 


mean-field from Eq. 
kn) where S is the dirac delta function, Eq. 


A.6 


S' 

reduces to: 


(E(r,z)) = lim {N/S) K^aesei 

S —V C / ^2! 


(A.8) 


The reflection coefficient of the coherently reflected light from the layer of particles is thus given 
by: 


nTE (E(r,zo)) ^ 1 

R = -FT? -r = -Oeff £1 I — 1 - 

Ei(r,zo) 2 \cJ kzi 


(A.9) 


where (E(r, zq)), and Ei{r, zq) are the complex scalar magnitudes of the vector quantities (E(r, zq)), 

N 

Ei(r,zo) respectively and p = lims_,.oo -w is the number of particles per unit area. To find the 

o 

coherently transmitted beam below the layer of particles, thickness of which is less than the skin 
depth of the material, we take the field to be continuous across the layer of particles m i-e-, 

f™ = (l-HR™) (A. 10) 

where is the transmission coefficient of the incident TE wave. The unknown factor K can 
be determined by considering the field surrounding the layer of particles. Since we have a field 
(1 -I- .R"'"®)Ei(r, z) surrounding the particles this shoul d als o be the illuminating field in Eq. 


i.e., we get (1 -I- IP^) = K, which, on inserting in Eq. 


A.9 


A.7 


gives us: 




^ /wy £1 Oeff 

^2 V c / kzi 


1 - P; 


'Loy £lQeff 
^ C / kzi 


(A.ll) 
























and inserting Eq. |A.ll|in Eq. |A.10|we get: 


1 


^TE 



O' 


kzi 


(A.12) 


Eqs. A.ll| and |A.12| are equivalent to the coefficients obtained using the effective index model in 
Sec. m 


Appendix B. Fresnel reflection and transmission coefficients for TM polarization 


Here, we give the form of the Fresnel reflection and transmission coefficients to be used in Eqs. 16 


and 17 when a TM polarized wave is incident on a layer of particles which is modeled as a thin-film 
using Maxwell-Garnett theory as described in Sec. Since the dielectric property of the film, as 
given from Eq. 15 has di ffere nt c omp onents along x and z directions the expressions for r, 
and to be used in Eqs. 16 andare given by pO] : 


kzlj^1 “f 


^(2) 


2kzi/ei ^ I £i 

kzlOl “b ^eff/^efT,a: y 

kzlj^l “b ^eff/^efT,a: y 


(B.l) 

(B.2) 

(B.3) 


where 


^eff 




^efF,3: 

^eff,z 


(B.4) 
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